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Dynamics of the normal-superconductor phase transition and the puzzle of the 

Meissner effect 
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The analisis of Pippard [1] for the growth of the normal phase into the superconducting phase 
in the presence of a magnetic field H > He is applied in reverse to the case H < He {He =critical 
magnetic field). We carry out the analysis both for a planar and a cylindrical geometry. As the su¬ 
perconducting phase grows into the normal phase, a supercurrent is generated at the superconductor- 
normal phase boundary that flows in direction opposite to the Faraday electric field resulting from 
the moving phase boundary. This supercurrent motion is in direction opposite to what is dictated 
by the Lorentz force on the current carriers, and in addition requires that mechanical momentum 
of opposite sign be tranferred to the system as a whole to ensure momentum conservation. In the 
cylindrical geometry case, a macroscopic torque of unknown origin acts on the body as a whole 
as the magnetic field is expelled. We argue that the conventional BCS-London theory of super¬ 
conductivity cannot explain these facts, and that as a consequence the Meissner effect remains 
unexplained within the conventional theory of superconductivity. We propose that the Meissner 
effect can only be understood by assuming that there is motion of charge in direction perpendicular 
to the normal-superconductor phase boundary and point out that the unconventional theory of hole 
superconductivity describes this physics. 

PACS numbers: 


I. INTRODUCTION 

When a metal is cooled into the superconducting state 
in the presence of a magnetic held H, the magnetic held is 
expelled from the interior. A current hows within a Lon¬ 
don penetration depth (Al) of the surface that generates 
a magnetic held that exactly compensates the magnetic 
held in the interior. This is the Meissner effect, discov¬ 
ered experimentally in 1933 [2]. In this paper we analyze 
the dynamics of this process and argue that it cannot be 
understood within the conventional BCS-London theory 



FIG. 1: Three possible routes for the magnetic field ex¬ 
pulsion in a cylindrical superconductor. The dots represent 
magnetic field lines coming out of the paper. The arrows give 
the direction of the currents (7). See text for a discussion of 
the processes. 


of superconductivity [3, 4], contrary to what is generally 
believed. 

First let us discuss some general features of the pro¬ 
cess. Figure 1 shows three conceivable routes to get from 
the initial to the hnal state of a long cylinder cooled into 
the superconducting state in the presence of a uniform 
magnetic field along the cylinder axis. For simplicity we 
assume translational invariance in the direction of the 
cylinder axis (long cylinder). In route (a), a current de¬ 
velops near the surface of the cylinder that gradually and 
uniformly depresses the magnetic field in the interior, un¬ 
til the magnetic field vanishes. In terms of magnetic field 
lines, their density decreases uniformly in the interior 
of the cylinder. In (b), small superconducting domains 
nucleate at random positions in the sample (determined 
by imperfections or temperature fluctuations), each with 
its own surface current that cancels the magnetic held 
in its interior. These domains grow and coalesce with 
each other, when domains coalesce their interior currents 
cancel out and only the boundary currents remain. Even¬ 
tually one single domain that occupies the entire sample 
results. In route (c), a single domain grows from the 
center, with a surface current that nullihes the magnetic 
held in its interior, and its boundary expands radially 
outward until it reaches the boundary of the cylinder. 

Note that in this cylindrical geometry the magnitude 
of the currents around the domains in the scenarios (b) 
and (c) is always the same, assuming it hows within a 
London penetration depth of the surface of the domain, 
since it always nullihes the magnetic held in the interior. 

We argue that the scenario (a) is untenable on physical 
grounds. In this scenario, while the transition is taking 
place a non-zero magnetic held exists throughout the in¬ 
terior preventing the establishment of phase coherence in 
the interior region, hence preventing formation of the su- 
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FIG. 2: (a) Growth of the normal (n) phase into the super¬ 
conducting (s) phase for a cylinder, under application of a 
magnetic field H larger than the critical field H^- (b) The 
reverse process (Meissner effect), under an applied magnetic 
field H smaller than the critical field. For both cases, the 
magnetic held at the normal-superconductor boundary is the 
critical held He- The shading of grey indicates the magnitude 
of magnetic held, white color indicates no magnetic held. 



FIG. 3: Growth of the normal (n) into the superconducting (s) 
phase for a planar geometry, under application of a magnetic 
held H = + p) , with p > 0. The magnetic held points 

out of the paper, the shading of grey indicates schematically 
its intensity. The n-s phase boundary (dashed line) moves 
down. 


ever these treatments do not address the fundamental 
questions of interest here. We will return to this point in 
a later section of this paper. The treatment in this paper 
applies to type I superconductors. 


perconducting phase that would lower the system’s free 
energy. In the absence of condensation energy there is 
no source of energy to drive the surface current and pay 
the electromagnetic energy price involved in changing the 
magnetic field in the interior. Thus scenario (a) can be 
discarded on purely theoretical grounds. 

Experimentally, there is evidence that the transition 
can occur through the process described in (b) [5] and 
also in (c), if a slightly smaller magnetic field is imposed 
at the center of the sample [6] . We argue that there is no 
fundamental difference in the processes (b) and (c) and 
that to understand the conceptual issues of interest here 
it is sufficient to concentrate on the scenario (c) which is 
what we will do for the rest of this paper. 

In a seminal paper in 1950, Pippard [I] analyzed the 
process of growth of the normal into the superconducting 
phase when a magnetic field larger than a critical field is 
applied, as shown schematically in Fig. 2(a). He argued 
persuasively that the rate of propagation of the phase 
boundary is governed by electromagnetic processes, and 
calculated the speed at which the boundary propagates, 
in approximate agreement with experiment. The applied 
field is H > He and it is reduced to He at the normal- 
superconductor boundary through the current generated 
due to Faraday’s law. Here we will apply the same rea¬ 
soning to the reverse process shown in Fig. 2(b), i.e. the 
Meissner effect. We will find that understanding the dy¬ 
namics of this process requires making assumptions that 
are incompatible with the conventional BCS-London the¬ 
ory of superconductivity. For completeness, we repeat 
Pippard’s analysis of the process (a) in the next section, 
before proceeding to the case of interest (b). 

There have been more recent theoretical studies of the 
growth of the superconducting phase into the normal 
phase [7, 8], taking into account the possibility that the 
planar interface may become unstable and more compli¬ 
cated growth patterns develop, and studying the differ¬ 
ences between type I and type H superconductors. How- 


II. GROWTH OF THE NORMAL INTO THE 
SUPERCONDUCTING PHASE 


Instead of a cylindrical geometry we will consider the 
planar geometry shown in Fig. 3, following Pippard [1]. 
The physics is the same and the calculation is consider¬ 
ably simpler. Later in sect. IX we return to the more 
interesting case of a cylindrical geometry. 

The applied magnetic field (in the z direction) at the 
boundary of the sample {x = 0) is given by 

H{x = 0) = Heil+p) (la) 


with p > 0. As the phase boundary moves in the negative 
X direction, eddy currents are generated that generate a 
magnetic field opposite to the applied one that reduces its 
magnitude, thus limiting the speed at which the bound¬ 
ary moves into the superconducting phase. The phase 
boundary is at the point x = xo(t) where the magnetic 
field is exactly the critical field He- 


Hixo) = He 


(lb) 


In the normal phase (ccq < cc < 0), the equations to be 
considered are 

^ 47r ^ dH 47r 

V X H = —J ==> = - Jy (2a) 

c ax c 


- - IdH dEy 1 dH 


J = aE ==> Jy = (jEy 


(2c) 


Here, H = Hz, J = Jyij, E = Eyjj. Application of 
Faraday’s law using a contour with one edge at the phase 











3 


boundary and the other deep into the superconducting 
phase where H{x) =0 yields 


Ey[x = xo) = 


He dxo 


(3a) 


and replacing Ey by Eq. (2c) and using Eq. (2a) the 
condition Eq. (3a) is 


dH ina dxo 

-^h=^o = - — 


(3b) 


Taking the x-derivative of Eq. (2a) and replacing the 
right-hand side by Eqs. (2b), (2c) leads to 


d^H Ana dH 


dx^ 


dt 


We define 


y = 


Xo 


H = He{l + f{y)) 

and from Eq. (4), / satisfies the equation 
d^f ^df 


dy^ 


c2 ^°dt' 


( 4 ) 

(5a) 

(5b) 

( 6 ) 


Expressing the time derivative of / in Eq. (4) in terms 
of the time derivative of xq yields the equation 


f'iy) 

f'iy) 


2TTa d{xo{t)Y 


( 7 ) 


from which we conclude that the right-hand side is inde¬ 
pendent of time, hence Xq increases linearly with time. 
We write 


ape 




( 8 ) 


and the parameter a is determined by the boundary con¬ 
ditions, which are 


fiy = 0)=P 
f{y = l) = 0 

df _ 27rcr9(xo(t))^ 

dy’^='^~ c2 dt ■ 

Eq. (9c) follows from Eq. (3b). 

Replacing Eq. (8) into Eq. (7) yields 

/" -f apyf = 0 

f'iy) = -ape^^^~y^) 


(9a) 

(9b) 

(9c) 

( 10 ) 

(11a) 



FIG. 4: Parameter a determined by the condition Eq. (12) 
for growth of the normal into the superconducting phase, and 
for growth of the superconducting into the normal phase de¬ 
termined by the condition Eq. (29) of Sect. 111. The dashed 
lines are the values resulting from the linear approximation 
discussed in Sect. IX. 


rv 

fiy)=P-ap dy'e 
Jo 




and a is determined by the condition 
a 


[ dy'e^(^-y'^) = l. 

Jo 


(lib) 


( 12 ) 


For small p, a ~ 1 and it decreases as p increases, as 
shown in Fig. 4. Note that Eq. (8) implies that the 
higher the conductivity of the normal metal the slower 
the motion of the phase boundary, as one would expect. 

The current and electric held in the normal phase are 
hence given by 




Jy — . H^e 2 
47r Xo 


c dt 


dTTCT Xo C dt 


(13a) 

(13b) 


and the magnetic Held is given by Eqs. (5b) and (11b). 
i7, Jy and Ey all decrease in going from the surface of 
the material {y = 0) to the phase boundary [y = 1). 

In the superconducting phase Eq. (2a) is of course 
satisHed which we reproduce here for convenience. 


^ - 47r ^ dH Att 

V X iJ = J = = > =-J, 

c dx c 


V 


and the London equation [3] is satisHed 

dJ„ 


V X J= - 


47rA2 


H==> 


dx 


AttXI 


H 


(14a) 


(14b) 


which upon combining leads to 


dx^ 


Ai 


hence 


(15a) 
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d^Jy _ 1 

dx^ "Ai " 

hence 


H{x) = 


(15b) 


(16a) 



Jy{x) = Jyixg )e(^ (16b) 

The current at the phase boundary on the superconduct¬ 
ing side is, from Eq. (14a) 

~ 47rAL 

while on the normal side of the phase boundary it is given 
by, from Eq. (13a) 


( 2 ^ 0 ") 


cap 

A 

47r Xo 


(18) 


so that it is discontinuous. The electric field in the su¬ 
perconducting side is obtained from Faraday’s law Eq. 
(2b), the expression for the magnetic field Eq. (16a) and 
the time derivative of xq Eq. (8) 


Ey{x) 


Jx-xo)/Xl 

c at 


(19) 


Note that the electric field is continuous at the phase 
boundary (cf Eq. (13b)), as expected. 

Upon taking the time derivative of Eq. (16b) and using 
Eqs. (17) and (19) we find 


dJy _ 

at ~ AttXI ^ 


( 20 ) 


which has a simple interpretation. With Jy = ngqvy, 
with Us the density of superconducting carriers of charge 
q, Vy their velocity, and using the standard expression for 
the London penetration depth [4] 


1 Aiinsq^ 

\\ mc^ 


( 21 ) 


with m the mass of the carrier, Eq. (20) is 


at ~m ^ 


( 22 ) 


which describes free acceleration of carriers of charge q 
and mass m due to the electric field generated by Fara¬ 
day’s law due to the moving phase boundary. 

The same analysis can be applied to a cylindrical ge¬ 
ometry, except that in that case an analytic solution is 
no longer possible. However an approximate treatment 
which is accurate for not too large values of p shows that 
the results are very similar to the results discussed in this 
section [1]. 

Thus, the physics of this process, shown qualitatively 
in Fig. 3 and described quantitatively by Eqs. (8)-(22),is 


FIG. 5: Growth of the superconducting (s) into the normal 
(n) phase for a planar geometry, under application of a mag¬ 
netic field H = Hc{l — p) with 0 < p < 1. The magnetic 
field points out of the paper, the shading of grey indicates 
schematically its intensity. The n-s phase boundary (dashed 
line) moves up. 


clear and intuitive. The moving phase boundary changes 
the magnetic flux and generates the electric field Ey both 
in the normal and the superconducting regions. In the 
normal region the electric field produces eddy currents, 
and in the superconducting region it produces a freely 
accelerating current. At the phase boundary the current 
drops in magnitude by a large amount in going from the 
superconducting to the normal side. This is of course 
because scattering processes set in, described by the re¬ 
sistivity p = 1/cr. When the boundary moves slightly 
down, the supercurrent at the boundary drops discontin- 
uously as scattering processes set in, and its momentum 
is transmitted to the solid as a whole. 

Instead, as we discuss in the next sections, there is no 
similarly simple physical interpretation for the reverse 
process, when the superconducting phase grows into the 
normal phase. 


III. GROWTH OF THE SUPERCONDUCTING 
INTO THE NORMAL PHASE 

We now apply the same analysis to the reverse pro¬ 
cess, where the applied magnetic field is lower than the 
critical field and hence the superconducting phase grows 
into the normal phase, as shown in Fig. 5. Again the 
kinetics is determined by electromagnetic processes. As 
the phase boundary moves into the normal phase, the 
current enhances the magnetic field at the phase bound¬ 
ary, which is again the critical field Etc- As before, the 
speed of motion of the phase boundary is determined by 
the eddy currents generated in the normal phase. 

The boundary conditions are now 

H{x = Q) = (23a) 

H{x = Xo) = He (23b) 

with p > 0. Eqs. (3), (4) and (5a) still hold, Eq. (5b) 
gets replaced by 

H = He{l - f{y)). (24) 

Eqs. (6) and (7) still hold, and Eq. (8) is replaced by 

xoit)^ =R^-^t. (25) 
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assuming the initial distance between the phase bound¬ 
ary and the surface of the material is R. Eqs. (9a) and 
(9b) still hold, and Eq. (9c) gets replaced by 

df 2'Kad[xQ{t)Y 

= 

Eqs. (10), (11) and (12) get replaced by 

/" - apyf = 0 (27) 


X Jy(X,t) Jy(X,t+At) Ey 



FIG. 6: Transition from the superconducting state to the 
normal state in a boundary layer. The charge carriers change 
their speed from a large Vs to a small Vn given by Eqs. (35). 


f'{y) = 1) (28a) 


f{y) =P - ap [ ^) (28b) 

Jo 

and a is determined by the condition 

a C dy'e^^y'^-^) = l. (29) 

Jo 

For small p, a ~ 1, now as p increases a also increases, 
as shown in Fig. 4. 

The current and electric held in the normal phase are 
given by 


(30a) 

“ 47r xo c dt 


J-'y — - -tJcC — <~x, ^ 

Attg Xq C Ot 


(30b) 


and the magnetic held is given by Eqs. (24) and (28b). 
H, Jy and Ey all increase in going from the surface of 
the material {y = 0) to the phase boundary {y = 1). 

In the superconducting phase Eqs. (14) and (17) still 
hold. In particular, the current on the superconducting 
side of the phase boundary is given by Eq. (17) which 
we reproduce here for convenience 


at the phase boundary. In contrast to the previous case 
where the electric and magnetic helds decrease continu¬ 
ously in going from the normal to the superconducting 
phase, here both the electric and magnetic helds achieve 
their maximum values at the phase boundary and de¬ 
crease from there both in going into the normal and into 
the superconducting phase. 

Note that Eqs. (20) and (22) still hold here, i.e. 


_ 

dt 

dvy 

dt 


- ^ E 

- 47TXI " 

(34a) 

= -Ey 
m 

(34b) 


however they have a different physical interpretation that 
in the previous case. Here, Eq. (34) describes the slow¬ 
ing down of the carriers in the superconducting phase 
by the Faraday held that applies a force in direction op¬ 
posite to the motion of the carriers. The origin of the 
motion of the charge carriers and resulting current in the 
superconducting phase is not described by the electro¬ 
magnetic processes discussed in this section, in contrast 
to the situation in the case of the growing normal phase. 
We discuss the situation in detail in the next sections. 


IV. ENERGY AND MOMENTUM 
CONSERVATION IN THE S N TRANSITION 


while on the normal side it is given by 

so the currents at the boundary now run in opposite di¬ 
rections, as shown schematically in Fig. 5. The electric 
held on the superconducting side is still given by Eq. 

(19), 

Ey{x) = (33) 

which now points in opposite direction to Eq. (19) be¬ 
cause the boundary is moving in opposite direction. The 
electric and magnetic helds are of course again continuous 


Figure 6 shows schematically a boundary layer that 
changes from superconducting to normal in a time inter¬ 
val At. The carriers in the layer of thickness S = xo{t)dt 
change their speed from a large Vs when they are in the 
superconducting phase to a small Vn when they become 
normal carriers, namely (from Eqs. (17) and (18)) 


c H, 
AttXl qus 


(35a) 


Vn 


= ap 


c H, 
47rxo qUs 


(35b) 


hence there is a change of kinetic energy and of momen¬ 
tum per carrier, given by 


A 1 2 1 2 

Aek = -mv„ - -mv. 


(36a) 
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Cooper pair 


V. 


X impurity 
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FIG. 7: Schematics of the transition shown in Fig. 6. A 
Cooper pair is moving with a high center of mass speed Va 
initially. One of its members gets attached to an impurity, 
and as the pair breaks up the original kinetic energy of the 
pair is used up to pay for the binding energy. When the pair 
unbinds essentially all the initial kinetic energy was used up 
in the energy required to unbind the pair. In the process, 
the momentum of the pair is transferred to the impurity and 
hence to the lattice as a whole. 


Apy = m{vn - Vs). (36b) 


How are these changes compensated? 

The kinetic energy of the supercurrent is not dissipated 
in collisions, since the transition is reversible [9]. The 
lowering of kinetic energy occurs because the carriers go 
from the superconducting to the normal phase, hence 
have to pay the price of the superconducting condensa¬ 
tion energy. The change in kinetic energy per unit volume 
is, neglecting the very small kinetic energy in the normal 
phase 

1 

AEk = UsAck =-ns-mv^ = (37) 

Z OTT 

where we have used Eq. (35a) for Vs and Eq. (21) for 
Al. This is precisely the condensation energy per unit 
volume. We can think of the condensation energy as 
the binding energy of Cooper pairs. To unbind the pairs 
requires to pay this energy cost, and it is paid by the 
kinetic energy of the carriers which consequently slow 
down to essentially zero kinetic energy. The process is 
shown schematically in Fig. 7. It involves interaction 
with some impurity or defect, and in the process the mo¬ 
mentum originally carried by the pair is transferred to 
the lattice as a whole, with no irreversible loss of energy. 

Note also that the origin of the carrier’s velocities Vg 
and Vn in Eq. (35) is well understood. Vg originates in 
the action of the electric field Ey on the superconducting 
carriers. For x deep into the superconducting region the 
carriers are initially at rest. As the phase boundary ap¬ 
proaches to a distance of order A l the electric field starts 
to grow and drives the carriers according to Eq. (22). 
Assuming at time t = 0 the phase boundary is at a:o = 0 
and at time tg it reaches the position x, i.e. xo(to) = x, 
we have for the carriers’ speed at x 



qHcXl 


(l_e^A^) (38) 


Jy(X ’0 A Jy(x,t+At) Ey 

◄- n I V ◄- 

- s 


FIG. 8: Transition from the normal state to the supercon¬ 
ducting state in a boundary layer. The charge carriers change 
their speed from a small Vn in the positive y direction (assum¬ 
ing g > 0) to a large Vs in the negative y direction, given by 
Eqs. (39). 



where we used Eq. (19) for the electric field. The term 
gx/Ai -g negligible for \x\/\l » 1 and Eq. (38) is the 
critical velocity Eq. (35a). On the normal side, the veloc¬ 
ity at the boundary is simply determined by the Faraday 
electric field and Ohm’s law. 

As we discuss in the next section, there is no com¬ 
parable understanding of the reverse process whereby a 
normal layer becomes superconducting. 


V. ENERGY AND MOMENTUM 
CONSERVATION IN THE N -> S TRANSITION 


In the reverse process, where a boundary layer changes 
from normal to superconducting, the direction of the cur¬ 
rent changes when the layer becomes superconducting, as 
calculated in Sect. Ill and shown in Fig. 8. Hence the 
direction of the motion of the carriers reverses, and in 
addition the speed of the carriers increases by a large 
amount, from Vn to Vg, given by 

Vn = —ap - - (39a) 

dTTCCo qUg 


c H, 
4ttXl qUg 


(39b) 


The change of kinetic energy and of momentum per car¬ 
rier are given by 


Acfc = imuf - ]^mvl (40a) 

Apy = m{vg - Vn). (40b) 

For the change in kinetic energy per unit volume we have 

1 

AEk = UgAek = Ug-mvl = —^ (41) 

Z oTT 

just like in the previous case {s ^ n transition) but 
with opposite sign. The increase in kinetic energy is sup¬ 
plied by the condensation energy. Again the process is 
reversible [9]. 

There are several questions that need to be answered: 


me 
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(1) What is the physical origin of the speed Vs, Eq. 
(39b)? 

(2) How is the momentum change given by Eq. (40b) 
compensated? 

(3) What is the physical mechanism by which the con¬ 
densation energy becomes kinetic energy of the supercur¬ 
rent, as given by Eq. (41)? 

We argue that the conventional BCS-London theory of 
superconductivity offers no answers to these questions. 
Let us discuss each of these questions. 


(1) What is the physical origin of Vs? 

The Faraday field Ey (Eq. (30b) or (33)) points in the 
positive y direction (see Fig. 8). Carriers in the normal 
state move in the positive y direction, driven by Ey (as¬ 
suming g > 0), and when the boundary layer becomes 
superconducting the carriers suddenly reverse their di¬ 
rection and move with large speed Vs in the negative y 
direction, opposing the force qEy that drives them to 
move in the positive y direction. 

How does the condensation process, or the Cooper pair 
formation process, make the carriers acquire a center of 
mass momentum in direction opposite to the applied elec¬ 
tric force? Furthermore, this center of mass momentum 
is only acquired if a magnetic field is present: if a ma¬ 
terial goes superconducting in the absence of a magnetic 
field, no Vg is acquired by the carriers condensing into 
the superconducting state. 

Perhaps it will be argued that in the process of be¬ 
coming part of the superconducting condensate and es¬ 
tablishing phase coherence with the already supercon¬ 
ducting carriers, the carriers in the boundary layer that 
is becoming superconducting are ‘carried along’ by the 
carriers in the layer right below, that are moving with 
speed Vg. However this argument has two problems: (1) 
it does not explain how the speed Vg is generated in the 
first place, when the transition to superconductivity be¬ 
gins; (2) in the cylindrical geometry of Fig. (2b), the 
momentum (angular momentum in this case) of the cur¬ 
rent near the s — n boundary increases as the boundary 
moves out because its perimeter increases, so a simple 
momentum transfer from one layer to the next cannot 
account for this growth of total momentum. 

Thus, we argue that this momentum (or angular mo¬ 
mentum) of the normal carriers becoming superconduct¬ 
ing has to originate in an external source, not in the 
superconducting carriers themselves. The only conceiv¬ 
able external source is the magnetic field E[. However, we 
can conceive of no mechanism for E[ to impart momen¬ 
tum in the y direction to the carriers within conventional 
BCS-London theory [10]. 


(2) How is momentum conserved? 

There is no momentum of the electromagnetic field in 
the y direction (since it is in direction E x H) and as¬ 
suming the process is sufficiently slow no momentum is 
carried away by electromagnetic waves, so the only way 
to compensate for the change in momentum Eq. (40b) is 
to transfer it to the lattice. If we assume for definiteness 
that the carriers have charge g > 0 so they move in the 
same direction as the current, the lattice as a whole has 
to acquire momentum in the positive y direction (to the 
left in Fig. 8) when the boundary layer in Fig. 8 becomes 
superconducting. 

It may be argued that there is a mechanism based on 
the electron-phonon interaction to account for such a mo¬ 
mentum transfer to the lattice. This has not been dis¬ 
cussed in the literature, and we can conceive of no way in 
which the argument could be made. Furthermore there 
are a large number of materials classes for which it is 
generally believed that the electron-phonon interaction 
is not involved in the mechanism of superconductivity 
[11], where superconductivity is believed to be caused 
by electron-electron interactions. How can carriers that 
bind into Cooper pairs through a non-phonon mecha¬ 
nism interact with the lattice in such a way that they 
transfer the right amount of momentum, given by Eq. 
(40b), whose magnitude is determined by the magnetic 
field, to the lattice? We can conceive of no answer to 
this question within the conventional understanding of 
superconductivity [12]. 


(3) How is the condensation energy converted into 
kinetic energy of the supercurrent? 

The energy of the carriers is lowered by condensa¬ 
tion into the superconducting state. This energy is 
converted into the kinetic energy of the carriers in the 
boundary layer becoming superconducting as given by 
Eq. (41). However the energy lowering, presumably as¬ 
sociated with Cooper pairing, has no directionality. How 
does it get converted into a kinetic energy that is asso¬ 
ciated with motion in one particular direction, i.e. the 
(—y) direction in Fig. 8? Again we argue that the con¬ 
ventional BCS-London theory offers no answer to this 
question. 

In the absence of answers to these three questions 
we argue that the Meissner effect remains unexplained 
within the conventional theory of superconductivity. 


VI. WHAT IS WRONG WITH THE 
CONVENTIONAL UNDERSTANDING 

The questions that we raise here have not been raised 
in the context of the conventional theory of superconduc¬ 
tivity. In this section we try to understand why. 


Let us start by recalling the argument frequently used 
to make London’s equation Eq. (14b) plausible [13, 14]. 
Starting from Newton’s equation for dissipationless mo¬ 
tion of a charged fluid under an electric field E and using 
Faraday’s law, 


dJ _ ^ 

dt m 


(42a) 


dt 


(V X J) 


Usq"^ OH 
me dt 


(42b) 


and integrating Eq. (42b) in time, Eq. (14b) results 
“outside of a constant of integration” [13], with Ai given 
by Eq. (21). However, integrating Eq. (42b) in time 
yields 


2 

Vx J(t)-Vx J(t = 0) = -^(7L(t)-i?(t = 0)) (43) 

me 

where the r dependence is implicit. If a normal metal is 
cooled into the superconducting state in the presence of a 
spatially uniform magnetic field iLo through its interior, 
the initial superfluid current J(t = 0) = 0 and Eq. (43) 
yields 

V X J(r, t) = - ^(i?(f, t) - TLo) (44) 
me 

which is not London’s equation Eq. (14b) . Quite the 
contrary, Eq. (44) implies that J{r,t) = 0 and ff(r, t) = 
Hq for all times t > 0, so that the magnetic field remains 
unchanged inside the superconductor. 

It is also argued that London’s equation and the Meiss¬ 
ner effect follow from BCS theory [4]. However this is 
only true to the extent that one assumes that in cooling 
a metal into the superconducting state the system will 
reach the state of minimum energy described by BCS the¬ 
ory, which implies the existence of a macroscopic wave- 
function 


^-(r=) = |^'(r)|e*^(”) (45) 

with a unique well-defined macroscopic phase 6{r) de¬ 
scribing all the electrons in the superfluid. From Eq. 
(45) it follows, with = |'I'(r)|^ and assuming uniform 
|4/(r)|, that the supercurrent is given by 

J=!M(W0--H) (46) 

m c 

hence Eq. (14b) results from taking the curl. However 
this is a circular argument. Assuming that Eq. (45) is 
valid implies that phase coherence has been established 
throughout the system, which implies that the magnetic 
field has been expelled, since no phase coherence can be 
established in the presence of a magnetic field. The ques¬ 
tion of how phase coherence is established, i.e. what are 
the dynamical processes that lead from the normal state 


with no phase coherence to the superconducting state 
described by Eq. (45) and the magnetic field expelled, 
needs to be addressed to make this argument valid. Oth¬ 
erwise, in the absence of experimental evidence, it could 
also be concluded that the BCS state described by Eq. 
(45) will never be attained by a system cooled below its 
critical temperature in the presence of a magnetic field. 

It is argued that time-dependent Ginzburg Landau 
theory (TDGL) [15, 16] describes the time evolution of 
the superconducting order parameter 41 and in particu¬ 
lar explains the Meissner effect, as assumed in the treat¬ 
ments of refs. [7, 8]. However, in TDGL it is assumed 
that the time evolution of the order parameter is such 
that it will evolve towards its equilibrium value that min¬ 
imizes the Ginzburg-Landau (GL) free energy. In partic¬ 
ular, that a generalized force exists, given by the func¬ 
tional derivative of the GL free energy with respect to 
the order parameter, that drives the time evolution of 
the order parameter, making it relax to its equilibrium 
value in roughly the time it takes a carrier to diffuse over 
a distance of a coherence length [17]. This is an assump¬ 
tion that cannot be rigurously derived from microscopic 
BCS theory, unlike the equilibrium GL free energy that 
can under certain approximations. In the ‘derivations’ 
of TDGL from BCS it is assumed that the system will 
reach the BCS ground state. For example, in Ref. [15] it 
is postulated that making the electron-electron attraction 
g time-dependent, i.e. g —>■ e^*g, i5 > 0, so that the inter¬ 
action is zero for t —^ —oo and attains its value g at t = 0, 
results in an order parameter A —e‘^*A. Again this as¬ 
sumes the answer, i.e. that the BCS state with expelled 
magnetic field will be attained starting from a normal 
metal in the presence of a magnetic held when the metal 
becomes superconducting, without proving that this is 
the case and without specifying the physical processes 
by which this occurs. 

More specihcally, within TDGL the current in the su¬ 
perconducting region is always given by Eq. (46), with 
ris = |\I'(r,t)|^. As a boundary layer becomes supercon¬ 
ducting the current grows not because the speed of the 
carriers changes continuously from to Vs but because 
the number of superconducting carriers Ug increases con¬ 
tinuously as |4>(r, t)| increases. The speed of the super¬ 
conducting carriers is hxed, at the value Vg = —{q/me)A 
determined by the London equation. Thus, in this de¬ 
scription the velocity of an individual carrier changes 
discontinuously in sign and magnitude from Vn to Vg 
when the carrier becomes part of the superconducting 
condensate described by 4'. The question of momentum 
conservation is not even addressed. We argue that this 
is not a satisfactory description of the physics. 


VII. WHAT IS REQUIRED TO EXPLAIN THE 
MEISSNER EFFECT 


We have argued in the introduction that the Meiss¬ 
ner effect necessarily involves motion of the normal- 
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X carrier becoming superconducting 

FIG. 9: As carriers become superconducting (crosses), they 
thrust forward in the x direction a distance \l and in the pro¬ 
cess acquire a speed Vg in the negative y direction (assuming 
positive carriers) due to the action of the Lorentz force. The 
normal carriers in the x-interval (zo, xq -|- An) undergo back- 
flow in the negative x direction due to electrostatics. In the 
process they acquire impulse in the positive y direction due 
to the Lorentz force and transmit this impulse to the lattice 
through normal scattering processes. 


superconducting phase boundary, namely that the sce¬ 
nario of Fig. 1(a) that does not involve motion of the 
phase boundary can be excluded on purely theoretical 
grounds. In order to answer the questions posed in Sect. 
V we propose that it is necessary to assume that the 
motion of the phase boundary is associated with motion 
of charge. In the conventional understanding of super¬ 
conductivity no motion of charge is associated with the 
motion of the phase boundary [7, 8]. 

Assume that as carriers condense into the supercon¬ 
ducting phase they experience a sudden thrust in direc¬ 
tion normal to the phase boundary and towards the nor¬ 
mal region, as shown schematically in Fig. 9. For sim¬ 
plicity we assume q > 0 for the moment. The Lorentz 
force acting on these carriers moving with velocity Vx in 
the positive x direction 


Fl = -Vx X H =--VxHcy (47) 

c c 

points in the negative y direction, and if it is larger than 
the Faraday force qEy it will impart motion to the car¬ 
riers in the negative y direction which is the direction of 
the current Jy in the superconducting region. The total 
force including the force from the Faraday field Ey at the 
phase boundary (Eq. (33)) is 

Fl = q{Eyy + ^ X H) = - ^{vx - ^)Hcy (48) 



FIG. 10: Schematic depiction of the orbits of the carriers in 
the theory of hole superconductivity. Normal state carriers 
have small non-overlapping orbits of radius supercon¬ 

ducting state carriers have large overlapping orbits of radius 
2Al. As normal carriers become superconducting their orbits 
expand, they acquire azimuthal velocity due to the mag¬ 
netic field, and charge is transferred from the superconducing 
to the normal side of the phase boundary. 


Therefore, the carriers becoming superconducting have 
to thrust forward a distance Xl into the normal region 
to acquire the required speed Eq. (39b). 

As these superconducting carriers are thrusting into 
the normal region in the positive x direction they will 
cause a counterflow of charge from the normal region in 
the negative x direction due to electrostatic forces, in 
order to restore charge uniformity, as shown in Fig. 9. 
The counterflowing carriers are still in the normal state, 
and experience scattering with lattice imperfections and 
phonons. Through the Lorentz force Eq. (47) they ac¬ 
quire the same impulse in the positive y direction (since 
Ux < 0 for these carriers) as the thrusting superconduct¬ 
ing carriers acquired in the negative y direction, and they 
transmit this impulse to the lattice as a whole through 
normal scattering processes. This accounts for the mo¬ 
mentum conservation question discussed in section V. 

Finally, the third question raised in section V about 
the conversion of condensation energy into kinetic energy 
of motion in a definite direction is also answered by this 
process. The energy for the thrusting forward is provided 
by the condensation energy, the direction of the current 
carrying the kinetic energy is determined by the Lorentz 
force Eq. (47) acting on the carriers thrusted in the +x 
direction. 


VIII. EXPLANATION OF THE MEISSNER 
EFFECT WITHIN THE THEORY OF HOLE 
SUPERCONDUCTIVITY 


so we can ignore the Faraday field livx » dxo/dt. The 
speed in the y direction that these carriers acquire if they 
thrust forward a distance Ax in a time interval At is 


Vy=r^dt = -^-^r\xdt = -^-^Ax (49) 

Jo rn c Jq c 

and for Ax = Al 


V 


y ~ 


-^XlH, = 
me 


c H, 
AttXl qus 


(50) 


The theory of hole superconductivity [18] proposes a 
concrete realization of the physical processes described 
in the previous section. Within this theory, charge carri¬ 
ers change their character from hole-like to electron-like 
when they pair and condense into the superconducting 
state [19]. In a semiclassical description, electronic or¬ 
bits expand from a microscopic radius {kp =Fermi 
wavevector) to mesoscopic radius 2Xl [20]. This orbit 
expansion is driven by lowering of kinetic energy [21], 
and has associated with it expulsion of negative charge 
from the interior to the surface of superconducting bodies 
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[22]. In the geometry considered here, it leads to trans¬ 
fer of negative charge from the superconducting into the 
normal region at the phase boundary, as shown schemat¬ 
ically in Fig. 10. 

That the system in the normal and superconducting 
states can be characterized by orbits of radius and 
2Xl respectively can be seen from the magnetic suscep¬ 
tibility. The Larmor diamagnetic susceptibility for elec¬ 
trons of density n per unit volume in orbits of radius r 
perpendicular to the magnetic field is 

XLarmor[r) = - -jC (51) 

This yields the Landau diamagnetic susceptibility of the 
normal state when the orbits have radius kp^: 

XLarmorij' = kp ) = ~ ( 52 ) 

with g^ep) = 3nl2ep the density of states and gp = 
|e|?i/2meC, and perfect diamagnetism when the orbits 
have radius 2Xp 

TIG^ 1 

XLavmorir = 2 Xl) =-^^^^{2Xl)‘^ =- — . (53) 

The fact that electrons occupying large orbits can ac¬ 
count for the perfect diamagnetism of superconductors 
was pointed out by Slater [23] and others long ago. 

As the orbit expands to radius 2Xl in the presence 
of a magnetic field, the carrier in the orbit acquires an 
azimuthal velocity [24] 


gAi 

V0 = - He 

me 


(54) 


due to the action of the Lorentz force. This is precisely 
the speed Eq. (39b) (using Eq. 21), and thus provides 
a dynamical explanation of the Meissner effect [20]. Be¬ 
cause the electrons have negative charge the azimuthal 
motion acquired is in the counterclockwise direction. The 
superposition of these orbits gives rise to net electronic 
charge motion in the +y direction, corresponding to cur¬ 
rent Jy in the negative y direction as given by Eq. (31). 
The rate of radial expansion is given by Vx defined in 
the previous section, which can be plausibly assumed to 
be much larger than dxo/dt. As the orbits at the phase 
boundary expand, they will enclose normal small orbits 
right above them (Fig. 10). These normal electrons will 
be pushed in the negative x direction due to electrostatic 
interactions, since the expanding orbits carry negative 
charge with them that extends into the region of the nor¬ 
mal phase right above the phase boundary. This accounts 
for the backflow discussed in the previous section. 

Thus, the theory of hole superconductivity provides a 
realization of the conditions discussed in the previous sec¬ 
tion that are needed to account for energy and momen¬ 
tum conservation as the phase boundary advances into 
the normal region in the presence of a magnetic field. 



FIG. 11: Superconducting region expanding from the center 
in a cylindrical geometry. 0 < r < ro is the superconducting 
region, ro < r < i? is the normal region. The magnetic field 
points out of the paper in the z direction. 


IX. CYLINDRICAL GEOMETRY 


Next we consider the growth of the superconducting 
phase in a cylindrical geometry as shown in Fig. 11. 
Here, the paradox associated with momentum conserva¬ 
tion (angular momentum in this case) is even more ap¬ 
parent than in the planar case discussed earlier. 

The magnetic field points along the axis of the cylinder 
{z direction), and the boundary between the supercon¬ 
ducting and normal phase is given by radius r = rQ{t), 
so that the system is superconducting for 0 < r < rg 
and normal for xq < r < R, with ro(t = 0) = 0 and 
T'oito) = R, with to the time it takes for the magnetic 
field to be expelled, to be determined. The induced Fara¬ 
day field as the superconducting phase expands from the 
center, Eg{r), points in the +9 direction, generating an 
azumuthal current Jg in the normal region that creates 
a magnetic field in the +z direction that tries to restore 
the magnetic field being expelled. The analogous of Eqs. 
(2) for the cylindrical geometry are 


^ ft 47r ^ dH 
y X H = —J ==> 

c or 


Att 

=- Jg (55a) 


c at r or 


—cUt 


J = aE ==> Jg = (jEg (55c) 

and the boundary condition analogous to Eq. (3a) is 

A.(co) = f^. (56) 
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Replacing Eg by Eq. (55c) and using Eq. (55a) the 
condition Eq. (56) is 


dH ^ dTTcr dvQ 

^ 

and combining Eqs. (55) yields 


(57) 


similar to Eq. 
conditions are 


1 d dH _ AiradH 
r dr dr dt 

(4) for the planar case. 


(58) 


The boundary 


H{r = R) = Hc{l-p) (59a) 



H{r = ro) = H, (59b) 


and Eq. (56). 

In the planar case, separation of variables was achieved 
by defining the variable y = x/xq, and the boundary 
conditions were given at y = 0 (x = 0) and y = 1 
{x = xq). This doesn’t work in the present case. If 
we define y = rjrg, the boundary conditions are given at 
y = I and y = R/xq, which are not time-independent. 
Time-independent boundary conditions can be achieved 
by dehning r = {xq — R)y + R, so that the conditions Eq. 
(59a) and (59b) are given at y = 0 and y = 1 respectively, 
but in terms of this variable Eq. (58) does not separate 
as Eq. (7) did. 

Eq. (58), with boundary conditions Eqs. (57), (59) can 
of course be solved numerically. Instead, we discuss here 
an approximate analytic solution, following Pippard[l], 
which becomes exact in the limit p —>■ 0. To zeroth order 
in p we may assume that the magnetic field in the normal 
region is given by 77 = 77^, which yields for the induced 
electric field in the normal region 


Eg{r) 


He rp drp 
c r dt 


(60) 


Using Eqs. (55c) and Eq. (55a) integrated from xq to R 
yields 


dxQ R 
ro 


dt 


9 

pc 

dTTCT 


(61) 


and performing the time integration from 0 to t with 
ro(t = 0) = 0 and ro(t) = Tq yields 


(!£)2[i + 2Zn-] = ^ (62a) 

it ro to 


with 


io 


iraR^ 

pc^ 


(62b) 


the total time for the magnetic field to be expelled from 
the cylinder. 


FIG. 12: Radius of the superconducting region ro{t) versus 
time t 


These results are exact in the limit p —?> 0. In the next 
order of approximation we assume a linear interpolation 
for the magnetic held between x = xq and x = R. This ap¬ 
proximation, for the planar case, yields the dashed lines 
shown in Fig. 4, corresponding to a = 3/(3-|-p) for the ‘n 
growing’ case and a = 3/(3 — p) for the ‘s growing’ case. 
It can be seen in Fig. 4 that that they closely follow the 
exact solutions for the planar case, hence we may sim¬ 
ilarly assume that the results will be close to the exact 
results in the cylindrical case. Hence we assume 


77(r) = 77,(1-p^^) (63) 

R- Xq 

varying linearly from 77, at r = ro to 77,(1 — p) at r = 77. 
This leads instead of Eq. (60) to 


Eg[x) — -tT ^ ^ ra - 12 

c X dt xo{R — xo) 


r!L _ lo + 

I 3 3^2^° 


(64) 


Using Eqs. (55a) and (55c) we hnd that the time required 
for held expulsion now is slightly less than given by Eq. 
(62b), 


to = -(65 

pc^ 9 

The equation giving the time evolution of xq in this ap¬ 
proximation is lengthy and will not be reproduced here. 
In Fig. 12 we plot the time evolution of xq for various val¬ 
ues of p obtained within this approximation, the curve la¬ 
beled p = 0 corresponds to Eq. (62a). It can be seen that 
the results obtained for the time evolution for different 
p’s are very similar provided the time is scaled according 
to the relation Eq. (65). From this we conclude that the 
exact results obtained from numerical integration of Eq. 
(58) would be very similar to the results predicted by the 
lowest order approximation Eq. (62). 
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These results are valid in the normal phase. In the 
lowest order approximation, the current at the n-s phase 
boundary on the normal side is given (from Eqs. (60) 
and (55c) by 


Je{ro) 


cp 

Airroln^ 


He 


( 66 ) 


and flows in the positive 9 direction (counterclockwise). 

In the superconducting phase the fields and current are 
governed by Ampere’s law and London’s equation [3] 

47r ^ dH 

V X /L = — J ==> = - Js (67a) 

c or c 


V X J = - 


dTTAi 




4n\l 


H. (67b) 


Combining Eqs. (67a) and (67b) yields 


i_a ^ _ J_ 


( 68 ) 


The electric field in the superconducting region is ob¬ 
tained from Faraday’s law Eq. (55b) and the expression 
for the magnetic field Eq. (71) 


1_9 

r dr 


{rEe{r)) 


H{r) drp _ ^ 
cA L dt 2ro 


(74) 


with the variation of ro{t) with t given by Eq. (61). As¬ 
suming the dominant variation comes from the exponen¬ 
tial behavior in Eq. (71), which is the case for r >> Ai, 
yields for the electric field in the superconducting region 


Ee{r) = 


ro)/\L jT^dro \l 

cXl \ r dt 2r 0 


so that for Al << ro,r 


dJeir) 

dt 



Eg{r) 


(75) 


(76) 


holds as in the case of the planar interface, describing the 
decelaration of the carriers in the superconducting region 
due to the electric field that points in direction opposite 
to the superconducting current flow. 


and the solution with boundary condition 77 (rg) = He is 
[25, 26] 


Angular momentum 


H = He 


loir/Xr) 

loiro/XL) 


(69) 


with Io{x) a Bessel function of imaginary argument, with 
asymptotic behavior 


Io{x) ->• —7== (70) 

y l-Kx 

for large x. Hence, assuming rg >> Xr the magnetic 
field in the region r < rg is given by 


As the superconducting region expands, the current 
carriers in the superconducting region acquire increas¬ 
ing angular momentum, both because the angular mo¬ 
mentum of each carrier of the supercurrent near rg in¬ 
creases proportionally to rg and because the number of 
carriers within A^ of the phase boundary grows as the 
perimeter of the circle of radius rg increases. It is found 
experimentally[27-29] that the mechanical angular mo¬ 
mentum carried by the Meissner current for a cylinder 
of radius R and height h with applied magnetic field H 
parallel to the cylinder axis has magnitude 


i7(r) = (71) 

and from Eq. (67a) the current in the superconducting 
region is given by 

and flows in the negative 9 direction (clockwise). Hence 
at the n-s phase boundary on the superconducting side 
the current is (for rg >> Al) 

* = (73) 

i.e. it is much larger in magnitude than the current in 
the normal side (Eq. (66)) and flows in opposite direc¬ 
tion. Thus, just like in the case of the planar interface 
there is a large discontinuity in magnitude and a change 
in sign of the current at the phase boundary, as shown 
schematically in Fig. 11. 


Le = ^R^hH (77) 

2|e| 

and points in direction parallel to the applied magnetic 
field. In Eq (77) nie is the bare electron mass. This 
relation holds for all superconductors [27-29] and results 
if the carriers in the Meissner current are bare electrons 
of mass TOg and negative charge e [30]. 

We can ignore the angular momentum carried by the 
eddy currents generated in the normal phase, because it 
is much smaller than the angular momentum carried by 
the superconducting carriers. When the phase boundary 
is at radius rg the total electronic angular momentum is 
then 

Reit) = -^^hHero{tfz (78) 

2e 

pointing in the positive z direction. By conservation of 
angular momentum we have to assume that the body as a 
whole acquires an equal and opposite angular momentum 

Lbodyit) = ^^^hHeroit^z (79) 
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and given the time evolution of tq, Eq. (62a), we have 

1 


Uodvit) = 

27re tol + 21n— 


(80) 


with V the volume of the body. This corresponds to an 
angular velocity of the body 


}{t) = 


rripC Hr t 


7r|e|p tg 1 -(- 2ln^ 


(81) 


with p the body’s mass density. Thus, as the supercon¬ 
ducting phase grows, a macroscopic angular momentum 
and angular velocity of the body as a whole is generated 
that grows approximately linearly with time, given by 
Eqs. (80) and (81). The torque acting on the body that 
causes its rotation is, from Eq. (79) 


dL, 


body 


dt 


dro ^ 

- hHcVo—z 

e dt 


(82) 


which we can write using Eqs. (61) and (62b) as 


nieC He hR^ „ 

e 4to In— 
ro 


(83) 


For example, for a sample of i? = 1cm, h = 5cm, Hr = 
200G, if it takes to = Is to expel the magnetic field, 
T = 1.4 X grem^/s^ when tq = R/e = 0.37R. 

The conventional theory of superconductivity provides 
no explanation for the origin of this macroscopic torque, 
nor for how the electronic angular momentum and ionic 
counter-angular momentum are generated and grow as 
the superconducting phase grows. In contrast, the the¬ 
ory of hole superconductivity provides an explanation ac¬ 
cording to the physics discussed in Sect. VIII, shown 
schematically in Fig. 13. In a cylinder of radius rg and 
height h, the total angular momentum due to carriers of 
density Us moving in orbits of radius 2 Al with orbital 
speed vg is 


Le = [meVe{2XL)]ns{TTrl)h (84a) 

which can be rewritten as 


Le = [meVero]ns{2nroXL)h. (84b) 


In the first form, it describes the aggregate angular mo¬ 
menta of carriers of density rig occupying an area tttq, 
each in an orbit of radius 2 Al with angular momentum 
meVe{2XL)- In the second form, it describes the angu¬ 
lar momentum from carriers in a ring of radius rg and 
thickness Al, i.e. the carriers of the Meissner current 
when the superconducting region has radius rg. Clearly 
these two descriptions are equivalent, since superposing 
the 2Xl orbits the internal velocities cancel out and only 
the surface current remains. The speed of the carriers 
at radius rg is given by Eq. (73). Assuming the carri¬ 
ers are bare electrons of density rig, as implied by the 
experiments [27-29], Eq. (73) with Jg = engVg yields 


Vg = - 


c 

dTrAirigC 


He 


(85) 




FIG. 13: Left panel: pattern of currents and fields in cylin¬ 
drical geometry. The magnetic field points out of the paper. 
As the superconducting phase grows, the boundary supercur¬ 
rent acquires increasing angular momentum. Assuming the 
current is carried by negative charge carriers (electrons), as 
shown by experiment [27-29], the body as a whole rotates in 
the clockwise direction to conserve angular momentum. The 
right panel depicts the expanded orbits in the superconduct¬ 
ing region that explain the process within the theory of hole 
superconductivity. 


and from Eq. (82) 


rUpC n 

Le = -^rlhHe (86) 

in agreeement with Eq. (78). 

The generation of angular momentum occurs as fol¬ 
lows: as the phase boundary moves out, the orbits of 
normal carriers becoming superconducting carriers ex¬ 
pand from microscopic radius kp^ to radius 2Xl- In the 
presence of magnetic field He, the electron in the expand¬ 
ing orbit experiences a Lorentz force and acquires veloc¬ 
ity Vg in the positive 9 direction given by Eq. (85), or 
equivalently Eq (54) [24], and hence angular momentum 
in the +z direction. The sum of these dynamically gener¬ 
ated angular momenta yields the total electronic angular 
momentum Eq. (82). Furthermore, as the orbits at the 
phase boundary expand, there is an outflow of negative 
charge from the superconducting into the normal region 
that gives rise to a radial backflow of normal carriers as 
discussed in Sect. VII. The carriers flowing inward ac¬ 
quire angular momentum in the —z direction through the 
action of the Lorentz force and transmit this angular mo¬ 
mentum to the lattice by collisions. In this way, as the 
phase boundary moves out, electrons acquire increasing 
angular momentum in the +z direction and the body ac¬ 
quires the compensating angular momentum in the —z 
direction given by Eqs. (78) and (79). 

We argue that there is no way to explain the dynam¬ 
ics of this process of angular momentum generation in 
the absence of the radial outflow and inflow of charge 
discussed here. 








14 


X. SUMMARY AND DISCUSSION 

In this paper we have calculated the electromagnetic 
fields and currents in the process where the normal phase 
grows into the superconducting phase and in the reverse 
process where the superconducting phase grows into the 
normal phase, in the presence of a magnetic field, as¬ 
suming for simplicity a planar interface, and also for a 
cylindrical geometry. We have also obtained the speed of 
motion of the phase boundary in these processes. 

The results obtained are easily understood for the pro¬ 
cess where the normal phase grows. However, for the re¬ 
verse process, which correspond to the Meissner effect, we 
have argued that the results cannot be understood within 
the conventional understanding of superconductivity, i.e. 
the London equation, BCS-Eliashberg theory and time- 
dependent Ginzburg-Landau theory. This conventional 
framework describes pairs of carriers that change their 
center of mass velocity discontinuously in both sign and 
magnitude as they bind into a Cooper pair and condense 
into the superconducting state in the presence of a mag¬ 
netic field, without providing a physical explanation for 
how this change of velocity, in direction opposite to that 
dictated by the force exerted by the Faraday electric field, 
occurs, and without explaining how the change of mo¬ 
mentum is compensated so that the physical law of mo¬ 
mentum conservation is not violated. It is expected that 
the change in momentum is compensated by momentum 
transferred to the lattice as a whole, but no mechanism 
for such transfer is provided and it is not clear how this 
would occur, particularly for superconductivity mecha¬ 
nisms that do not involve the electron-phonon interac¬ 
tion. 

For a cylindrical geometry, the question of momentum 
conservation is even more puzzling than in the planar 
geometry, because the momentum grows with time rather 
than staying constant in total magnitude as in the planar 
geometry. We calculated the rate of growth of angular 
momentum of the conducting carriers and of the body 
as a whole, and the magnitude of the angular velocity 
of the body and of the torque acting on the body as 
the magnetic field is expelled. These are macroscopic 
observables for which we argue there is no explanation 
within the conventional theory of superconductivity. 

It is true that for a microscopic quantum process, such 
as the transition of an electron from a 2p to a Is atomic 
orbital with emission of a photon, we do not object to 
the fact that the electron changes its velocity discontin¬ 
uously. However, the superconductor is a macroscopic 
system, albeit quantum, and in the spirit of Bohr’s cor¬ 
respondence principle one would expect that its behav¬ 
ior should be understandable without invoking processes 
that seem to fly in the face of macroscopic physical laws, 
as the TDGL description does. And even for microscopic 
atomic processes there is an account of energy and mo¬ 
mentum (or angular momentum) conservation (angular 


momentum is carried away by the photon in the above 
example), while no explanation of how momentum or an¬ 
gular momentum conservation is satisfied in the Meissner 
effect is provided by the TDGL formalism. Finally, the 
fact that one state of a system has a lower energy than 
another state does not predict that the system will evolve 
from the high energy to the low energy state, even in the 
microscopic realm. For example, two Be atoms will not 
spontaneously fuse to form an O atom. 

Instead, we have argued that the Meissner effect can 
be understood by making assumptions that are not part 
of the conventional understanding of superconductivity, 
namely that there is motion of charge in direction nor¬ 
mal to the phase boundary associated with the motion 
of the phase boundary. More specifically, that there is 
motion of “superconducting charge”, that undergoes no 
scattering, from the superconducting into the normal re¬ 
gion, and backflow of “normal charge”, that does undergo 
scattering, in the opposite direction, as the phase bound¬ 
ary advances into the normal region. To understand the 
Meissner effect it is not necessary to assume a specific 
sign for the charge involved in the charge flow perpendic¬ 
ular to the phase boundary. However within the theory 
of hole superconductivity the sign is uniquely defined: it 
is negative charge that moves from the superconducting 
into the normal region [31]. 

Note that the proposed processes resemble processes 
known to occur in superfluid '^He. The process where 
a superconducting region expands into a normal region 
will occur if the superconducting region is supercooled 
relative to the normal region. In superfiuid ^iLe, it is 
well known that there is flow of superfluid from colder to 
warmer regions and counterflow of normal fluid (fountain 
effect), precisely what is required to explain the Meiss¬ 
ner effect as discussed here. We have proposed that this 
commonality of behaviours occurs because both for su¬ 
perfiuid ^He and for superconductors described by the 
theory of hole superconductivity, the transition into the 
superfluid or superconducting state is driven by lowering 
of kinetic energy [32], rather than lowering of potential 
energy as predicted by the conventional theory. 

The Meissner effect has been known for 82 years, and 
it is generally believed that it is explained within the con¬ 
ventional BGS-London theory of superconductivity. We 
have argued here that this is not the case, and that the 
Meissner effect can only be understood if there is charge 
motion in the direction of motion of the phase boundary, 
which is not predicted by the conventional theory but is 
predicted by the alternative theory of hole superconduc¬ 
tivity. Thus we argue that the Meissner effect remains 
unexplained within the conventional theory, and hence 
that the conventional theory in its present form can only 
apply to superconductors that do not exhibit the Meiss¬ 
ner effect, in other words that it does not describe any of 
the currently known superconducting materials [11]. 
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